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1. Introduction 

The differential equations in the squared external momentum p 2 for the 
Master Integrals (Mis) of the 2-loop sunrise graph with arbitrary masses 
mi,?^,?^ of Fig. ^ were written in pQ. They were used for obtaining 




Fig. 1. The 2-loop sunrise graph. 



analytically particular values and behaviours at zero and infinite momentum 
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transfer pQ, at pseudo-thresholds {2] and threshold [3], as well as for direct 
numerical integration 4 . In this contribution I will report on some progress 
in the anlytic study of the solutions of the equations for arbitrary momentum 
transfer in the equal mass limit mi = 1. A more complete account will be 
given elsewhere (Sj; while the algebraic burden in the arbitrary mass case 
will surely be much heavier, there are indications 6 that the approach can 
be extended to the arbitrary mass case as well. 

In the equal mass limit the 2-loop sunrise has two Mis, which in the 
usual d-continuous regularization scheme can be written as 
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Let us put p 2 = z (z is positive when p is Euclidean); the two Mis then 
satisfy the following linear system of first order differential equations in z 
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The system can be rewritten as a second order differential equation for 
S(d, z) only 
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As the second MI S±(d, z) can be written in terms of S(d,z) and its first 
derivative 
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we can take from now on S(d, z) and its derivative dS(d, z)/dz as the effec- 
tive Mis. 
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2. From near 4 to near 2 dimensions. 

We want to expand S(d, z) around d = 4 as Laurent series in [d — 4) and 
then to obtain analytically the values of the coefficients of the expansion by 
solving the relevant differential equations. It was found a posteriori that all 
the formuale are much simpler when expanding around d = 2. To give the 
relations between the two expansions, let us recall that acting on any scalar 
Feynman integral in d dimensions with a suitable differential operator, one 
obtains the same integral in (d — 2) dimensions, times a numerical factor 
depending on d [Z|. Acting on the Mis in d dimensions (or, in our case, 
on the two functions S(d,z) and dS(d, z)/dz), one obtains the same Mis in 
d — 2 dimensions in terms of mass derivatives of the Mis in d dimensions, 
which can be expressed again in terms of Mis in d dimensions; solving the 
linear system for the ci-dimensional Mis and replacing finally d by d + 2 one 
obtains 
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Quite in general, if 



A{2 + d) = B{d) , 
one can set d = 2 + ry and Laurent-expand in 77; one obtains 

k k 

The Laurent expansion in r\ of S(d, z) for d = 4 + r] begins with a double 
pole in r\ and reads 

5(4 + v ,z) = (4, z) + -S^ (4, z) + (4, z) + rjS^ (4, z) + 

while S{2 + 77, z) has no singularities in 77, and its expansion is 

S{2 + r],z) = S(°)(2, z) + 7/5^(2, z) + 

By inserting the two expansions in Eq.©, one gets the required coefficients 
S^(4, z) of the Laurent expansion in 7/ of S(4 + r/,z) around 4 in terms 
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of the coefficients 5^(2, z) of the expansion of 5(2 + r],z) around 2. As 
S(d, z) is regular at d = 2, the poles in 77 are not hidden in S(d, z) but are 
explicitly exhibited by the l/(d— 2) factors in the r.h.s. of Eq.l|5j). Working 
out the algebra, one finds at once 
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3. The expansion at d = 2 of the differential equation. 

By expanding systematically in (d— 2) all the terms appearing in Eq.©, 
one obtains a set of chained equations of the form 
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where the homogeneous part is the same for any order k, and the first few 
inhomogeneous terms are 
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The equations Eq.lfS]) are chained, in the sense that the inhomogeneous term 
of order k involves lower terms, of order (k — 1) (for k > 0) and {k — 2) (for 
k > 1) in the expansion of 5(2, z), as can be seen from Eq.© and is shown 
explicitly in Eq.s©. 

The system Eq.(J7J) is to be solved bottom up in k, starting from k = 
(in which case the inhomogeneous term is completely known) and then 
proceeding to higher values increasing k by one, so that at each step the 
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inhomogeneous term is known from the solution of the previous equations. 
The chained equations can then be solved by using Euler's method of the 
variation of the constants. The homogeneous equation is the same for all 
the values of k, 
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if ^i(z), ^2 (-z) are two independent solutions of the homogeneous equation, 
W(z) the corresponding Wronskian 



W(z) = *i(z) 
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dz dz 

according to Euler's method the solutions of Eq.sQ are given by the integral 
representations 
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(k) (k) 

where ^\ ,^ 2 are two integration constants. 

Eg. 1)11(1 at this moment is just a formal representation of the solutions 
for the coeffcients S <yk \2,z); it becomes a substancial (not just formal!) 
formula only when all the ingredients - the two solutions of the homogeneous 
equation ^fi(z), their Wronskian W(z) and the two integration constants 

(k) 

are known explicitly. 
Although the Wronskian is defined in terms of the ^i(z), it can be 
immediately obtained (up to a multiplicative constant) from Eq.@. An 
elementary calculation using the definition Eq, (|l(J|) and the value of the sec- 
ond derivatives of the *$>i(z), as given by Eq.Q of which they are solutions, 
leads to the equation 

' + -^—) W(z) 



±W{z) 
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z + z + 1 
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which gives at once 
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where the multiplicative constant has been fixed anticipating later results. 
Finding the two ^i(z) requires much more work. 
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4. Solving the homogeneous equation at the singular points. 

By inspection, the singular points of Eq.@ are found to be 

z = 0, —1, —9, oo ; 

at each of those points one has two independent solutions, the first regular 
and the second with a logarithmic singularity. The expansions of the so- 
lutions around each of the singular points is immediately provided by the 
differential equation itself. 

Around z = the two solutions of Eq.@ can be written as 

*?>(*) = 

^\z) = lnz4 0) (z)+4 0) (z) , (13) 
where the ip^(z) are power series in z. Imposing ?/^(0) = 1, one finds 

4%) = -j^ + fW...; (14) 

the coefficients are given recursively (hence up to any reuired order) by the 
equation. The radius of convergence is 1 (the next singularity is at z = — 1) 
and the two solutions are real for positive z (spacelike momentum transfer). 
The continuation to the timelike region is done by giving to z the value 
z = —(u + ie); for < u < 1 one has Inz = Inu — iir and ^ij (z) develops 

an imaginary part — zV'i/'i (#)• 

Similarly, around z = — 1 the 2 independent solutions can be written as 

^\z) = ^\z) 

^\z) = ]n(z + l)^\z)+^\z) , (15) 

with Eq.(jUJ) providing recursively the coefficients of the expansions in powers 

of (z + 1) of the two ij)^{z) once the initial condition is given. If 4 X) (0) = 1 
one has 

$\*) = l + ^+l) + ^ + l) 2 + ... , 

4\z) = +|(* + l) + i|(* + l)2 + ... > (16) 

with radius of convergence 1 (up to the singularity at z = 0) etc. 

The other two singular points z = —9 and z = oo can be treated in the 
same way, the corresponding formualae are not given due to lack of space. 
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5. The interpolating solutions. 

Having the solutions piecewise is not sufficient, one must build two solu- 
tions in the whole — oo < z < oo range by suitably joining the above expres- 
sions of the solutions at the singular points. A hint is provided by the knowl- 
edge of the imaginary part of the original Feynman integral S(d,p 2 ) Eq.(^Q) 
at d = 2 dimensions; as already observed in [S], the Cutkosky-Veltman rule 
gives for the imaginary part of S(d,p 2 ) at d = 2 and u = — z > 9 (and up 
to a multiplicative constant) the integral representation 

/•(\/S-i) 2 db 
J ^ = / /B T m ' W 

J 4 y/R 4 {U,b) 

where R 4 (u,b) stands for the polynomial (of 4th order in b and 2nd order 
in u) 

R 4 (u,b)=b(b-4)(b-(V^-l) 2 )(b-(V^+l) 2 ) , 

and the b integration runs between two adjacent zeros of R 4 (u,b). As the 
inhomogeneous part of Eq.(j3J) cannot develop an imaginary part, the imag- 
inary part of the Feynman integral in d = 2 dimensions, J(u) of Ea.l)17[). is 
necessarily a solution of the associated homogeneous equation at d = 2, i.e. 
of Eq.Q - a fact which can also be checked explicitly. 

One is then naturally lead to consider all the 6-integrals of 1/ ^R 4 (u, b) 
between any two adjacent roots for all possible values of u. The details of 
the analysis cannot be reported here again for lack of space. As a result, 
one finds for instance that when u is in the range < u < 1 the roots of 
R 4 (u, b) are ordered as 

< (a/u — l) 2 < (v^+ l) 2 < 4 

and the associated ^-integrals are 

"(V^- 1 ) 2 db 



o y/-R 4 (u,b) ' 
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(0I-1) 2 y/R±(u,b) 

4 db 



(VS+i) 2 y/-R 4 {u,b) 



The three integrals are all real (and positive) due to the choice of the sign 
in front of R 4 (u,b) in the square roots; more important, they all satisfy 
Eq.@ - therefore they cannot be all independent. With standard changes 



of variables, they can be brought in the form of Legendre's complete elliptic 
integrals of the first kind 8 ; in that way one finds for instance 

3 ^/(1 + V«) 3 (3- V5) \0- + VQ 3 (3-y/u)) 

showing in particular that the first and third integrals are indeed equal. A 
similar (but different!) formula holds for the second integral. 

Although usually one can not do very much for expressing elliptic inte- 
grals in terms of other more familiar functions (such as logarithms or the 
like), in the limiting cases u —* and u — * 1 two of the 4 roots of R^{b,u) 
become equal, and the by now elementary 6-integrations give 



(18) 



(19) 

One has now all the information needed for defining two solutions ^i(z) 
of Eq.@ in an interval which contains the two singular points z = and 
z = — 1. Let us start by defining, for z > 0, 

= ^\z) , 
* 2 (*) = *< o) (a0 . 

That fixes the multiplicative constant in the Wronskian as well, giving the 
result already anticipated in Ea.(|12j). From Ea.s (|13ll4|) we easily read the 
behaviours of the ^i(z) for u = —z small and positive; but in the range 

< u < 1 the solutions can also be expressed in terms of the 4° ( u )l by 
matching the behaviours u — > + of the j| (u), Ea.s(|18j) to the behaviours 
of the ^i(z), one finds, in the interval < u < 1, 
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One can now compare Ea.s lfT9*)l with Ea.s H15ll6|) and express the 4°' 1} (, u ) 
in terms of the SH^^z)) substituting in Ea. (|2Uj) one finds for the solutions 
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^fi(z), for z around —1, the values 

4tt Air 
m 2 ( z -ie) = ^1. ^ln21n3-2vr-i91n2^S 1) (z-ie) 

+ ^(-21n3 + m)^ 1) (z-ie) . 

One can then move to the next interval 1 < u < 9 and so on till the ^%(z) are 

expressed, in the whole range — oo < z < oo, in terms of the ^[ k \z), each 
known within the convergence radius of the expansions given by Eq. Q , as 

well as in terms of the elliptic integrals J- k ' l \u). 

6. The integration constants. 

S(d,p 2 ) Eq.Q is known to be real for u = —p 2 below the threshold 
at u = 9. Take the solution as given by Euler's formula Ea. (|ll|) ; in the 
region < u < 1, or > z > - 1 the argument w of the inhomogeneous 
term N^ k \w) varies in the interval > w > z > —1, and is therefore 
real (N^ k \w), Eq.(jHJ) involves either real algebraic fractions or lower order 
terms of the expansion in (d — 2) of S(d, z), which are real in that region) - 
therefore an imaginary part, if any, can come only from the ^i(z) and the 
^/i(w). By using the values of the ^i(z) as given by Eq.s(J20J), one finds for 
u = —z in the range < u < 1 

lmS^(2,z) = -V3^ 2 k) 4 ' 1] (u) 

implying , for any k, 

The argument can be repeated in the interval 1 < u < 9 (between 
pseudothreshold and threshold), where the ^i(z), in analogy with Eq.s(J20J), 
are expressed in terms of 



J^\u) 
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One finds for z in the interval —l>z> —9, i.e. 1 < u < 9 

ImS< fc >(2, z) = S^J^i-z) ( *f ) + 2V3 f ' J^-J^\- W )N^ (2, 

TT \ JO W(W) 
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as the imaginary part must vanish in that interval, the other integration 
constant is given by 

7. Conclusions. 

The values of all the quantities entering in Ea. QllJI . namely the two 

*f?i(z), W{z) and the two integration constants i£j have been obtained, so 
that the previously formal expression given by Eq. became a substancial 
formula giving the functions S^ k \2, z) in closed analytic form. Indeed, from 
the explicit knowledge of the singularities and the relevant expansions of 
the *&i(z), the singularities and the relevant expansions of the S( fc )(2, z) are 
immediately obtained - and from the explicit knowledge of the singularities 
and the relevant expansions of the S^(2,z) the fast and precise numerical 
routines for their evaluation can in turn be obtained. 
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